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ABSTRACT 

In this article, we test the hypothesis that Cepheids have infrared excesses due to mass loss. We fit 
a model using the mass-loss rate and the stellar radius as free parameters to optical observations from 
the OGLE-III survey and infrared observations from the 2MASS and SAGE data sets. The sample 
of Cepheids have predicted minimum mass-loss rates ranging from zero to 10~*Mq yr~^, where the 
rates depend on the chosen dust properties. We use the predicted radii to compute the Period- 
Radius relation for LMC Cepheids, and to estimate the uncertainty caused by the presence of infrared 
excess for determining angular diameters with the infrared surface brightness technique. Finally, 
we calculate the linear and non-linear Period-Luminosity (P-L) relations for the LMC Cepheids at 
VIJHK + IRAC wavelengths and we find that the P-L relations are consistent with being non-linear 
at infrared wavelengths, contrary to previous results. 
Subject headings: Cepheids - circumstellar matter - Magellanic Clouds - stars: mass loss 



1. INTRODUCTION 

Classical Cepheids are powerful standard candles be- 
cause they follo w the Period-L uminosity (P-L) relation 
or Leavitt Law (jLeavittl I1908D . Th e P-L relation has 
been used as a tool for Galactic (iFeast &: WhitelodJ 
Il997| ). extragalactic (jGieren et al.l I2009L and references 
there i n) and cosm ological studies (iFreedma n et al.l 
[2OOIL iSandage et al . 1998). The LMC Cepheids have 
been central for deriving these relatio n s, initially us- 
ing o p tical data (lUdalski et al.l Il999al: ISandage et al.l 
200l IKanbur et all 120031 IKanbur fc Neeowl 120061: 
Fouaue et al.l I2007D . then near-infrared observations 
(iGroenewegenI 120001 : iPersson et all 120041 : iFouque et all 
I2007D. and recently using infrared observations 
(jFreedman eTHI 120081 : iNgeow fc Kmlbiul l2008h . The 
infrared P-L relations were determined by matching 
sources from the SAGE (Surveyi ng the Agents of a 
Galaxy's Evolution) archival data (jMeixner et al.ll2006l) 
with the known Ceph eids. 

Ngeow et al.l ()2009[ ) computed the infrared P-L rela- 
tions for a larger s ample of 1848 C epheids from the 
OGLE-III Database (jSoszvnski et al.li2008.) by matching 
those sources with the two published epochs of SAGE 
data, and also averaging the two epochs of infrared fluxes 
to bring the fluxes closer to the mean flux. They tested 
the P-L relations for non-linearity, and found that the 
IR P-L relations are linear, while the relations at wave- 
lengths shorter than iiT-band are non-linear. However, 
there are a number of Cepheids with infrared fluxes 
that are > 3 standard deviations brighter than predicted 
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by the P-L relations. Similar results we re found using 
AKARI N-band data (|Ngeow et al.ll2010[ ). It is impor- 
tant to understand the nature of this excess infrared 
brightness because the intrinsic infrared flux of Cepheids 
is less mctallicity dependent than optical P-L relations 
(jFreedman et al.l |2008() , which makes the IR P-L rela- 
tions powerful tools for extragalactic studies. The James 
Webb Space Telescope will be able to observe Cepheids 
in distant galaxies, making it possible to use the IR P-L 
relations to determine distances, but the unknown source 
of infrared excess will increase the uncertainty of the re- 
sults. 

One hypothesis for the source of the infrared ex- 
cess is the existence of a circumstellar envelope (CSE) 
of dust around the LMC Cepheids, an alogous to the 
CSEs observed around Galactic Cepheids (iKervell a et al.l 
,2006; Merand et al. 2006, 2007). Muhi- wavelength ob- 
servations of RS Puppis and / Car con firm the pres- 
ence of CSEs around these Cepheids (jKervella et al.l 
120091 ). There is more evidence of infrared ex- 
cess in Ga lactic Cepheids usi ng IR AS observations 
(jMcAlarv fc Welch 1986; Dcasv 11988f ) and Spitzer ob- 
servations ( Mareng o et al.l [20090 ■ On the other hand, 
Marengo et al. (2010l ) argued that there is no evidence 
from IR excess from warm dust based on Spitzer observa- 
tions and instead was due to carbon monoxide emission. 
The Cepheid RS Pup is also asso ciated with a nebula 
(lHavlenh972VKervena et al.ll2008D. 

KcrvcUa et al. (2006) speculated that CSEs surround- 
ing Galactic Cepheids were caused by mass loss tha t can 
leads to an infrared excess. iNeilson fc Lesteii ()2008f ) pro- 
posed a driving mechanism for this wind that is a combi- 
nation of radiative acceleration, pulsation and shocks in 
the atmosphere of the Cepheid. This proposed driving 
mechanism was tested by developing an analytic mass- 
loss model that was used to calculate rates of mass loss 
for a sample of Galactic Cepheids, with values between 
10~^° to \Q~'^ MQ/yr being found. This model predicts 
the fluxes of CSEs aro und the Cepheids t hat a gree well 
with those observed. INeilson fc Lesten (j2009( ) applied 
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this mass-loss model to theoretical models of Galactic, 
LMC, and Small Magellanic Cloud Cepheids and found 
that Cepheids in all three galaxies have significant mass- 
loss rates and infrared excesses. There is other evi- 
denc e for Cepheid mass loss in addition to infrared ex- 
cess; iNardetto et al. (2008 ') observed Ha line profiles of 
Galactic Cepheids with the purpose of detecting a hydro- 
gen circumstellar medium. They found that the Ha line 
profiles are asymmetric and blue-shifted in long period 
{P > 10 day) Cepheids, possibly caused by stellar wind. 

While it is important to characterize infrared excess 
of Cepheids in order to use them for various applica- 
tions, it is just as important to understand how mass 
loss affects the structure and evolution of Cepheids. 
Mass loss is a potential solution to the mass discrepancy 
of Cepheids, which is the difference between Cepheid 
masses estimated using stellar evolution calcula tions 
and estimates using stellar pulsation calculations (jCoxl 
|1980[ ). Currently, stellar pulsation calculations predict 
mas ses about 10-20% sma l ler than stell ar evolution mod- 
els IliliES^oi [200l EiieB [IOQI) . Measurements 
of dynamic masses of Cepheids in binary systems tend 
to ag ree with pulsation calculations (Evans et al. 2006, 
120081 ). suggesting the discrepancy is in the stellar evo- 
lution modeling. If Cepheids lose an average of 10~* to 
lO~^M0/yr, then mass loss could solve the mass discrep- 
ancy. 

Mass loss is hypothesized to create optically thin 
circumstellar dust shells that cause infrared excess in 
Cepheids but do not affect the visual extinction of the 
Cepheids. This excess should be seen in in frared observa- 
tions of the LMC from the SAGE survey (jMeixner et al.l 
l2006f ) for which there are currently two epochs of observa- 
tions. In fact, infrared excess has been d etected in a num- 
ber of evolved AGB stars in the LMC (|Viih et al.l 120091: 
iGroenewegen et aI]|2009D . In our earlier paper testing 
mass loss in LMC Cepheids, we found that mass loss 
may be important, with about 10% of OGLE-II Cepheids 
showing statistically significant evidence for dust shells. 
However, the uncertainty of mean infrared flux of the 
Cepheids is a limiting factor. The predicted gas mass- 
loss rates of the LMC Cepheids in that sample range 
from 10^^^ to lQ~^M0/yr assuming a dust-to-gas ratio 
of 1/250. 

Mass-loss rates in this range affect both the zero-point 
and the slope of the infrared P-L relations. When the 
infr ared excess w a s rem oved from the OGLE-II sam- 
ple, iNeilson et al.l (|2009f ) found that the scatter of the 
data was reduced. We also found that the 3.6 to 5.8 
lim relations are consistent with being non-linear while 
the 8.0 fim relation are linear, differing from the ob- 
served relations (Ngeow & Kanbur 2008). The infrared 
excess caused by mass loss appeared to have a greater ef- 
fect at shorter periods than at longer periods. However, 
the infrared stellar luminosity is an increasing function 
of period, making the infrared excess is less noticeable 
for long-period Cepheids, even for larger mass-loss rates 
(« 10-8 MQ/yr). 

These results suggest that the near-infrared fluxes, 
JHK, are potentially affected by mass loss. This 
would mean distance estimates using the Infrared Sur- 
face Brightness (IRSB) technique have uncertainty due 
to mass loss and infrared excess. The infrared surface 



brightness technique uses an empirical fit of the depen- 
dence of surface brightness on {V — K)q to determine 
the mean val ue and amplitu de of a Cepheid's angular 
diameter (Gier en et al.lll999l ). This angular diameter in- 
formation is then converted to a radius and a distance 
by using radial velocity observations to compute the am- 
plitude of the radius variation. This technique has been 
shown to be robust, but infrared excess from mass loss 
suggests that the color (V — K)o may be overestimated. 
The objective of this paper is to analyze fundamental- 
mode Cepheids in the LMC that have stellar fluxes from 
the OGLE-HI survey and infrared data from the 2MASS 
and from two epochs of the SAGE surveys, and test if the 
Cepheids have infrared excess. We determine mass-loss 
rates from infrared excesses determined from a correla- 
tion of these data sets by assuming a mass-loss model. In 
the next section we describe the method for testing the 
model. In section 3 we fit the model assuming zero mass 
loss, and the results of fitting the mass-loss model are 
given in section 4. In section 5, we test how the removal 
of the 5.8 and 8.0 ^m data affect the mass-loss model as 
well as how the unknown pulsation phase of the IR fluxes 
affects the predicted mass-loss rates. The statistical com- 
parison of the two models and the uncertainty of the fits 
due to the unknown pulsation amplitudes of the IR fluxes 
is given in section 5. Using the results of the models, we 
calculate the Period- Radius relation (section 6), and test 
the uncertainty of the IRSB technique due to mass loss 
(section 7). In Section 8 we compute linear IR P-L rela- 
tions from the data and compare the pr e dictions with the 
observ ed re sults of lNgeo\y et al.1 (|2009f ). iFreedman et al.l 
(|2008[ ). and lMadore et al.l (|2009^ in Section 9. In Section 
10, we test the P-L relations for non-linearity. 

2. METHOD FOR DETERMINING MASS-LOSS RATES AND 
RADII OF CEPHEIDS 

To model the mass-loss rates of LMC Cepheids, we em- 
ploy the sample of Cepheids from the OGLE-HI survey 
with V and /-band data. To obtain JHK and IRAC ob- 
servations, this sample is correlated with th e 2MASS and 
SAGE surveys. These data, from ' Ngeow et all |2009), 
have fluxes at up to nine wavelengths that have been 
corrected for extinction at all wavelengths. For Cepheids 
with measurements in both epoc h s of th e SAGE data, we 
follow the practice iNgeow et al.l ()2009() of averaging the 
IRAC fluxes to reduce the uncertainty due to pulsation 
phase and amplitude of brightness variation. 

Our operating assumption is that infrared excess is 
due to a dust shell formed at a large distance from the 
Cepheid that is optically thin at visible wavelengths, but 
other explanations are possible. For example, an infrared 
excess could be caused by blending of stars or by mis- 
taken associations. In our previous work, false associ- 
ations were possible because the average separation of 
SAGE and OGLE-II sources was about 0.77", with some 
Cepheids having separations as large as 1 to 2". Note 
that the uncertainty of the coordinate systems of the 
surveys are smaller than the separation: the average un- 
certainty of the O GLE-II coordinates are about 0.1 — 0.2" 
(|Szvmahskill2005f ). Using the OGLE-HI data, the mean 
separation with the SAGE sources is about 0.2" , and 
the mean separation with 2MASS sources is about 0.1". 
However, the calibration error of the OGLE-HI coordi- 
nates is only 0.06" (Udalski et al. 2008). The separation 



Mass Loss in LMC Cepheids II. 



1 1 1 1 1 

^^ ** * t * 
* * * ^ + 


2MASS - 


-*** *^*i*llf\:*f*** ", 


^fr4:*A^^***t , 


+-H- 

i-. * * 
**. * *+ 

+ ++ * 


1 


1 1 1 1 

SAGE - 

h \^^* *** 

msi^^*** 
t* * *\ 1 1 


^;*^#|^^ 
f / 


■■■ ■■ ,,-■ ''V> 


^1 



1 

o 

^ 0.8 

m 

c 0.6 

o 

2 0.4 

03 

§■ 0.2 


1 0.8 

r 0.6 
o 

2 0.4 

03 

S- 0.2 


0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

Log P(d) 

Fig. 1. — The separation of the sources from the 2MASS (top) 
and SAGE (bottom) surveys matched with the fundamental mode 
Cepheids from the OGLE-III survey. The dashed lines represent 
the separation that is the average plus 3cr deviation. 



of the sources is shown in Figure [H To guard against 
false associations, Cepheids in the sample with separa- 
tions greater than the average plus 3cr deviation have 
been removed. 

We model the mass-loss rate and radius of each 
Cepheid in t he sample using the same method as in 
iNeilson et al.l ([2009) with several improvements. The ob- 
served luminosity at each wavelength is assumed to be 
the sum of the stellar luminosity and the contribution 
due to the dust shell. Furthermore, we assume that the 
dust is composed of silicates instead of carbon because 
most Cepheids are typically not carbon rich. Assuming 
silicon instead of carbon decreases the predicted mass- 
loss rate by a factor of 2-3. The stellar luminosity is 



L^(Star)=47rR27rB^(Teff), 
and the contribution due to a dust shell is 



(1) 



L^ (Shell) = Btt 



< a2 > Md ^A 



Q. X 



< a-^ > pvd 

B,{Td)[l-W{r)]dr. 



(2) 



This luminosity is different by a factor of 47r^ from the 
luminosity used in the first article and is the correct 
derivation of the relation; this is an error in the first 
article and the resultant mass-loss rates are a factor 47r^ 
too large but the other results remain the same. In our 
first article, (Ncilson ct al. 2009), we checked how well 
a blackbody represents the Cepheid flux in the V- and 
/-bands. The difference is less than a few tenths of a 
magnitude. The dust shell luminosity is dependent on 
the ratio of the dust particle's surface area to volume, 
< a^ > / < a'^ >, the mean mass density of dust, p, 
the mass-loss rate, Md, and the velocity of dust, Vd- The 
variable Q^ is the absorption efficiency, which is treated 
more precisely in this work. Before, we assumed the ef- 
ficiency was « 2 and here we treat it as a function of 
wavelength, Q^ « 27r < a > /A, an approximate re- 



lation that is reasonable for A < lOfim. The dilution 
factor, W{r) is a measure of the radiative equilibrium 
of the gas and dust with stellar radiation. Previously, 
we assumed th at both the dust and g as had the same 
dilution factor. Ilvezic fc Elitzuri (|1997l ) argued that the 
temperature of the dust varies as ex r~^'^, but the tem- 
perature of the gas goes as ex r~^/^. In this wo rk, we 
use the representation of Ilvezic fc Elitzuil ()1997D . The 
dust is assumed to form at the condensation tempera- 
ture of 1200 K, which leads to a dust-condensation ra- 
dius, Rd « (i?,/2)(Teff/1200 K)"^. The dust particles 
are assumed to range in size from a — 0.005 to 0.25 
lim with a mean densi ty of 3.7 g/cm^, and to follow 
the iMathis et al.l ()1977[ ) distribution. The dust travels 
at a velocity similar to the escape velocity of order 100 
km/s. These assumptions are discussed in earlier works 
(jNeilson fc Lesteril200l INeilson et al.ll2009l ). 

The effective temperature of the Cepheids is given by 
the color (V — J) o, corrected both for extinction (see 
iNgeow et a l.' 2009, for details) and for the pulsation pe- 
riod from[Bcaulicu ct al. (2001). This leaves the Cepheid 
radius, i?*, and the dust mass- loss rate as free param- 
eters, and the dust mass-loss rate is converted to gas 
mass- loss rate by assuming a dust-to-gas ratio of 1 : 250, 
which is the canonical Galactic value scaled by the rela- 
tive metallicity of the LMC. All of the following figures 
and discussion regarding the mass-loss rates refer to the 
gas mass-loss rate unless specified otherwise. Therefore 
we have two free parameters to fit observed magnitudes 
at up to nine wavelengths for each Cepheid. We fit the 
model by calculating the minimum x^ value. This re- 
quires assuming both a distance to each Cepheid and an 
uncertainty of the fiux due to pulsation phase. We as- 
sume a distance modulus of 18.5 magnitudes to the LMC 
with an uncertainty of the distance being ±0.1 magni- 
tudes due to the thickness of the LMC (see INeilson et al.l 
120091 for details). A final improvement is to assign 
an uncertainty to the fl ux due to pulsation. Following 
Ngcow fc Kanbun (J2008), we set the uncertainty of the 
infrared flux at 1/3 the /-band amplitude, instead of as- 
suming a constant value of the uncertainty. 

The approach of the following two sections is flrst 
to establish a basis of comparison by determining the 
Cepheid radii assuming zero mass loss. Then we calcu- 
lated t he mass- loss rates of the Cepheids an d use the 
F-test (|Kanbur fc Ngeowll2004i INgeow et aLi r2009. and 
references therein) , to determine the quality of the flt for 
the mass loss and the zero mass loss cases, where the 
value of F is 



F ^ 



a?(iV-l)-ai(iV-2) 



(3) 



where trj and cr| are the uncertainties of the fit using the 
Cepheid radius only and then using both the Cepheid ra- 
dius and mass-loss rate, respectively, and N is the num- 
ber of data points flt. Therefore, we may calculate the 
mass-loss rates of the LMC Cepheids and test the valid- 
ity of the model. In the flt, we do not necessarily have 
information for all nine wavebands, therefore Cepheids 
are rejected if either the y or / is not present or if two 
or more of the IRAC bands are not present or if fewer 
than flve of the bands are present. Out of the 1848 fun- 
damental mode Cepheids in the OGLE-III survey, this 
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leaves 1398 Cepheids from SAGE epoch 1 observations, 
1387 from epoch 2 and 1552 Cepheids from the average 
of the two epochs. 

3. FITTING THE RADIUS TO THE DATA 

To establish a reference sample, we begin by assuming 
no mass loss and then use Equation [1] to fit the mean 
radius of each Cepheid. We fit the radius to each epoch 
of SAGE observations as well as the average of the fluxes 
for the two epochs. In Figure [U we show the value of 
X^ and the radius for each Cepheid in each sample. In 
the figure, those Cepheids with a large separation on the 
sky, as shown in Figure [H are removed; for epoch 1 this 
eliminates 62 out of 1398 Cepheids, for the second epoch 
107 of the 1387 Cepheids are removed, and for the aver- 
age of the two epochs 73 out of 1552 are removed. Only 
those Cepheids with a fit of x^ < 1-1 are shown in the 
plot of the predicted radius of each Cepheid. The x^ cut 
of 1.1 for radius-only fit is taken so that the slopes of 
the predicted V and /-band P-L relations agree within 
the uncertainties with the slo pes of the obs e rved V and 
/-band P-L relations given by iNgeow et al.l (|2009[ ). The 
connection between the structures of the P-L relations 
and the quality of fit of the modeled data is discussed in 
Section 9. Fitting only the radius to the observed fluxes 
appears to be a reasonable model. 

It is worth checking how much the predicted radius and 
X^ value changes for each epoch. In Figure|31 we show the 
relative difference of the x^ values (Left) and the relative 
difference between the predicted radii (Right). The rel- 
ative difference is defined as the difference between each 
epoch of observation divided by the predicted quantity 
from the average of the two epochs. The values of x^ are 
sensitive to the differences of the infrared flux because 
of the increased uncertainty of the infrared fluxes due to 
the unknown phase of pulsation at the time of the ob- 
servations and to an IR pulsation amplitude that can be 
up to 0.3 magnitudes. However, it is clear that the pre- 
dicted radius is not sensitive to the differences between 
the infrared fluxes at each epoch of SAGE observation. 
The maximum difference of radius ranges from about 3% 
at 35 /?0 (log P = 0.6) to about 2% for a Cepheid with 
/?, = 120/?Q (logP = 1.4) Cepheid. The radius depends 
on the total brightness and not the shape of the spectral 
energy distribution. Therefore, the radius is not sensitive 
to changes in the infrared flux due to pulsation. 

4. FITTING THE DATA WITH A TWO-PARAMETER 
MODEL 

We now expand the model by allowing the mass-loss 
rates to be greater than zero, now using both the Cepheid 
radius and mass-loss rate in the fits. In Figure |4j we 
plot the x^ computed for the two-parameter model, along 
with the predicted radius of each Cepheid with a value 
of x^ < 1.25. There are the same number of Cepheids 
in this sample with x^ < 1-25 as there are in the fit of 
the radius only. Compared to Figure [2l the x^-plots of 
Figure |4] for the two-parameter model are very similar 
but appear to have slightly more scatter in the fits of the 
long-period Cepheids. 

In Figure [5] we show the predicted mass-loss rates of 
the Cepheids with x^ < 1.25, predicted for each epoch of 
SAGE data. The predicted gas mass-loss rates range 
from zero to 10~*Mo/yr, but it must be noted that 



these are the m.inimum mass-loss rates. The percentage 
of Cepheids with non-zero mass-loss rates are 45% for 
the first SAGE epoch, 45% for the second SAGE epoch, 
and 48% for the average of the two epochs. For many of 
the Cepheids with zero mass loss, a small increase of x^ 
from the minimum would produce mass-loss rates of the 
order 10~^° MQ/yr. The result appears to suggest that 
the unknown pulsation phase causes the infrared excess, 
which is discussed in the next section. The uncertainty 
of the mass loss spans orders of magnitude. The gas 
mass-loss rate depends on a number of assumptions, in- 
cluding the dust velocity, the dust-to-gas ratio, and the 
size and shape of dust particles. The assumed values for 
the dust-to-gas ratio, the dust velocity, and the dust size 
are all likely to be lower limits. If the properties of any of 
these three are different from the assumed values then the 
mass-loss rates may be significantly larger. For instance, 
we assumed that the dust-to-gas ratio is 1 : 250, which is 
the m aximum possible value; we noted in iNeilson et al.l 
()2009() that the ratio may be about 1 : 1000 or even 
smaller if the LMC metallicity is not simply scaled down 
from the Galactic metallicity. This suggests that the gas 
mass-loss rates may be larger by the same ratio. Also, 
the assumed distribution of dust grain sizes has an effect 
on the predicted mass-loss rates. The dust grain sizes as- 
sumed in this work are based on the model for Galactic 
Cepheids (jNeilson fc Lesteiil2008[ ). where the metallicity 
is larger. In the LMC, with a smaller metallicity, it is 
likely the dust grain sizes have not grown as large. If so, 
then the corresponding mass-loss rates need to be larger 
to match the observations. 

Because the two SAGE epochs recorded the Cepheids 
at random pulsation phases, we computed the values of 
X^, radius and mass- loss rates change for both epochs 
to assess how these variables change with phase; Figure 
[6] shows the percentage differences of these variables as 
a function of pulsation period. The relative change of 
X^ ranges from 0.1% to 200%, suggesting that the fit is 
sensitive to the phase of pulsation, but this range is the 
same as found for the x^ values for fitting the radius only. 
The radius variation between the two epochs is typically 
only a few percent but some of the fits show a larger 
variation than in the radius-only model. The mass-loss 
rate shows the greatest variation, from about 1% up to 
10^%, approximately four orders of magnitude. There 
are some Cepheids that have a mass-loss rate variation 
<< 0.1% and some with a mass-loss rate variation that 
is >> 10"^%. The first situation is where the predicted 
mass-loss rates are identical in both epochs whereas in 
the second case the mass-loss rate is effectively zero in 
one of the epochs but is significant in the other epoch. 
The number of Cepheids in both cases is small. We con- 
clude that the unknown pulsation phase leads to a large 
uncertainty in the mass-loss rate because the change in 
infrared luminosity (in magnitudes) causes an exponen- 
tial change in the mass-loss rate. 

We use the F-test to compare the radius-only model 
and the radius-plus-mass-loss model. If the value of F 
is > 4 then we can state with 95% confidence that the 
mass-loss model provides a better fit than the radius only 
model. For the first SAGE epoch F = 2000, for the 
second SAGE epoch F = 2800 and for the average of the 
two epochs F = 3200. This strongly suggests that mass 
loss is an important phenomenon in the LMC Cepheids. 
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Fig. 2. — (Left) The x^ fit of the radius to each Cepheid in the sample for each epoch of SAGE observations plus the average of the two, 
with large separation Cepheids removed. The dashed line in each plot represents a value of x^ = l-l- (Right) The best-fit radius of each 
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Fig. 3. — (Left) The percentage difference of x^ from the radius-only model for the fits of each epoch of SAGE observation. (Right) The 
percentage difference of the best-fit radius of each Cepheid for each epoch of SAGE observation. 



5. DEPENDENCE OF MASS LOSS ON THE UNKNOWN 
PULSATION PHASE 

We have shown that slightly more than half of the 
Cepheids in the sample are consistent with zero mass loss, 
and this fraction is what one would expect if the infrared 
excess was entirely due to the unknown phase of the ob- 
serve IR fluxes. This expectation can be tested by vary- 
ing the IR flux of a sample of 200 Cepheids (the first 200 
in the sample, which is effectively a random period distri- 
bution). The amplitude of the IR flux for each Cepheid 
is about 1/3 the I-band amplitude from OGLE-III, and 
for each of the 200 Cepheids the IRAC fluxes are allowed 
to vary within the amplitude. We compute the best-fit 
mass-loss rate, and radius for each Cepheid with seven 
different phases and count how many of these phases 
predict significant mass-loss rates, M > 10^^"* MQ/yr. 



These seven phases create eight bins, the first being 
where a Cepheid has a mass-loss rate < lO"^** M^/yr at 
every phase, the second bin is 1 where a Cepheid is found 
to have significant mass loss for only one of the seven 
phases and so on to the eighth bin where a Cepheid has 
a mass-loss rate > lO"""^^ MQ/yr for all seven phases. 

We plot the distribution of the Cepheids with signifi- 
cant mass-loss rates (> 10~^^ MQ/yr) in Figure[7]as well 
as the average period of the Cepheids in each bin. The 
distribution of Cepheids appears to be highest in the last 
bin, with the fraction decreasing with decreasing bin un- 
til the bin representing zero mass loss at all phase with 
which suggests that Cepheids predicted to have signifi- 
cant mass-loss rates at one phase tend to have significant 
mass-loss rates at all phases. If the infrared excess was 
completely due to the unknown phase then one would 
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1 and 2 SAGE observations. 



expect that the distribution would have approximately 
equal peaks at each end of the distribution, i.e. for zero 
mass loss and mass loss at all phases. The distribution 
would appear parabolic, with a maximum at zero and 
the minimum at bins 3, and 4, with another maximum 
at bin 7. However, this behavior is not seen here. The 
fraction of Cepheids with mass loss at all phases is al- 
most three times more than the fraction with zero mass 
loss at all phases and the distribution is skewed towards 
more mass loss. This suggests that the mass-loss rates 



are real and not an artifact of the unknown pulsation 
phase, although the unknown pulsation phase is clearly 
an important uncertainty. 

The average period of the Cepheids with mass loss 
at all phases r aises an interesting question, where < 
logP >= 0.48. INgeow et al] ()2009h cut off the 5.8 /xm 
data at logP = 0.536 and cut off the 8.0 /xm data at 
logP = 0.702. The average period of the Cepheids in all 
bins except the bin is less than the 8 iim cut-off period 
while bin 7 has an average period less than the 5.8 iim 
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cut-off period. This is not surprising; iNgeow et al.l (|2OO90 
applied the cut-off periods to remove Cepheids near the 
hmiting magnitude of the observations at these wave- 
lengths. This result does suggest we test the propensity 
for mass loss with 5.8 and 8.0 /im data removed from the 
calculations. 

Our first test recalculates the best fit radius and mass- 
loss rates after removing the 8.0 /im data from the sam- 
ple, and compares this to the mass-loss rates calculated 
using fluxes from all wavelengths; the comparison is 
shown in Figure [H The data are plotted for two groups: 
those Cepheids with logP < 0.702 and those with longer 
periods. The mass-loss rates of both populations have 
similar scatter when the 8.0 /xm fluxes are excluded. 
From this we conclude that photometric errors in the 
8.0 /im fluxes at the detection limit do not significantly 
affect the mass-loss model. Short-period Cepheids with 
8.0 /im b rightness of about 1 2 to 11 magnitudes (see Fig- 
ure 3 of INgeow et al.l (|2009l )) have mass-loss rates that 
vary by about a factor of only two when that flux is 
used to constrain the model. It is these Cepheids that 
appear to have the most significant IR excesses, yet the 
mass-loss rates are not affected by photometric errors. It 
is possible that the short-period Cepheids with 8.0 ^m 
brightnesses of about 14 to 13 mag. are affected by pho- 
tometric errors at the limiting magnitude but there is 
similar scatter for longer period Cepheids. It is not pos- 
sible to quantify how much this effect contributes, if any, 
but we can state that the contribution does not affect 
the results significantly. 

Next we recalculate the best fit radius and mass-loss 
rates with both the 5.8 and the 8.0 /im data removed, 
and we have grouped the data into two bins, those with 
periods logP < 0.536 and those greater. Again, we com- 
pare to the radii and mass-loss rates calculated using all 
the wavelengths; the results are shown in the right panel 
of Figure [51 The scatter in the mass-loss rates is similar, 
and the exclusion of the 5.8 and 8.0 /im data has not 
made a substantial reduction in the number of Cepheids 
predicted to have significant mass loss rates. The appar- 
ent rows at low mass-loss rates in Figure [8] are due to the 
stopping the calculations when it is clear the the best-fit 
mass-loss rate in the analysis is less than 10~^^MQ/yr, 
which is equivalent to zero mass loss in the analysis. The 
effective temperature is determined by the period and 
(V — I)o, and the radius is most dependent on the opti- 
cal and near-infrared fiuxes, which implies that the stellar 
luminosity is well determined. The IRAC fiuxes provide 
the strongest constraint on the mass-loss rates, indicat- 
ing that the 3.6, 4.5, and 5.8 /im fluxes help offset this 
effect at 8.0 /im. When we compute the IR P-L relations, 
we will reconsider the influence of the 8.0 /im flux limit. 

6. PREDICTIONS OF THE PERIOD-RADIUS RELATION 

Because one of the free parameters of our model is the 
stellar radius, the fitting procedure determines the ra- 
dius of each Cepheid in the sample, enabling us to create 
a Period-Radius (P-R) relation for LMC Cepheids. We 
derive two P-R relations, one using the simpler model in 
which the radius as the only free parameter, and a sec- 
ond P-R relation in which both radius and mass loss are 
included. Table [1] gives the slopes and zero-points of the 
P-R relations for each of the two models and for each 
epoch of the SAGE data. The P-R relation is expressed 



TABLE 1 

Best Fit Parameters for Predicted Period-Radius 

Relations 



Epoch 


Slope Zero Point 


Dispersion 


N 


F 


Radius Model 


1 
2 
l-h2 


0.690 ±0.003 1.129 ±0.002 
0.688 ±0.003 1.131 ±0.002 
0.689 ±0.003 1.130 ±0.002 


0.017 
0.017 
0.017 


692 
700 
817 


3.57 
0.66 

0.88 


Radius + Mass-Loss Model 


1 
2 
1+2 


0.692 ±0.003 1.127 ±0.002 
0.688 ±0.003 1.130 ±0.002 
0.690 ±0.003 1.128 ±0.002 


0.017 
0.017 
0.017 


701 
710 
822 


3.16 
0.63 
0.74 



in the form log i?(i?0) = alogP(d) -I- b determined us- 
ing a_3a_J_teratiyecli2ping routine, similar to that used 
bv INgeow et al.l ()2009l ). For the fit using just the ra- 
dius, the slopes and zero-points of the P-R relations for 
various combinations of SAGE epochs all agree within 
the uncertainty. Similarly, the slopes and zero-points of 
each P-R relation using the mass-loss model for the three 
SAGE epochs agree. A comparison of the predicted P- 
R relations from the radius-only model to the P-R re- 
lations found using the mass-loss model shows that the 
slopes and zero-points agree within the uncertainty, al- 
beit the zero-point of the radius-only model is marginally 
larger but well within the uncertainty due to distance to 
the LMC. Note that the unknown phase of the infrared 
observations does not affect the fit of the P-R relation 
because the affect of phase is random for each Cepheid 
and cancels out in the fit. 

Table |2l gives t he P-R re lations from several other in- 
vestigations. Gieren et al.l ([1993) used the Infrared Sur- 
face Brightness technique to d e termi ne radii of LMC and 
SMC Cepheids, iGroenewegenI (|2007D derived a P-R rela- 
tion from five galactic Cephei ds with measured dist ances 
and angular diameters, while JKervella et al.l (I2004D used 
the su rface brightness method as well. iRuoppo et al.l 
(|2004f l used the CORS method (a modifi ed version of 
the B aade- Wesselink method derived by iCaccin et al.l 
(|1981h ) and lGieren et al.l (|1998[) used the IRSB technique 
to find logi? = 0.750(±0.02 4) logP-H.075(' ±0.007) that 
is identical to the results of lLanev fc Stobia (Il9 95') using 
the surface brightness method. (jBono et al.ll998M erived 
P-R relations from their theoretical calculations for 
two different mass-luminosity (M — L) relations; the 
shallower M — L relation the larger slope for the P- 
R relatio n. Our resu l ts ag r ee well with the rela tions 
found by'Gicrcnct al| (|1999( ). lRuoppo et al.l ()2004[ ). and 
Groencwegcn (2007[), but not with the relations with 
slopes > 0.750. Our results using both the one- and 
two-parameter fits predict Cepheii radii that are consis- 
tent with results from observations and with the P-R 
relations. 

As an experiment, we tested our sample for the pos- 
sibility of non-linear P-R relations, define as two linear 
relations with a period break at P — 10 days, similar 
to non- linear P-L relations defined in iKanbur fc Ngeowl 
(|2004bf) . There is no physical reason to expect that the 
relations should be non-linear, but it is a useful check 
of the fitting models. We do not list the non-linear P-R 
relations, but have tested the possibility using the F-test 
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Fig. 7. — (Left) The number of Cepheids that have significant mass loss (> 10" 
where the IR flux varies. (Right) The average period of the Cepheids in each bin. 



2 3 4 5 6 7 

Significant IVIass Loss Bins 
'"'^^ Mq/pt) for a given number of phases of pulsation, 





in 










Q 


10 






i 


10 


o 




CO 


in 










;=; 


10 






5 












S 


10 


cr: 




(0 


in-^ 










i» 


10-^ 


m 




> 






10- 




10 = 



Log P < 0.702 
Log P > 0.702 




D D DO-H 



I D □ □mB]-it=-tEi ' 



D H -B-IH- 



10-^' 10"^" 10"'" ... 10" 10"'" 10"' 
Mass-Loss Rates for All Bands 



10"' 
10"= 
10"° 
10"" 
10"" 

10"^^ 
10"='^ 
10"=' 
10"== 




B m cibp Dm ,B n 



Log P < 0.536 
Log P > 0.536 



10"'" 10"" ... 10" 10" 10 10 10" 
Mass-Loss Rates for All Bands 
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again. Because there are more than 1000 Cepheids in 
each sample, a value of F > 3 indicates the relations 
are consistent with being non-linear. The results from 
the SAGE epoch 2 and average data are consistent with 
a linear P-R relation but the P-R relations found using 
the mass-loss model and radius-only model from the first 
SAGE epoch data are consistent with being non-linear, 
with F = 3.16 for the mass-loss model and F — 3.57 
for the radius-only model. This anomaly is interesting 
but not sufficient evidence for non-linearity of the P-R 
relation. 



7. 



THE EFFECT OF MASS LOSS ON THE INFRARED 
SURFACE BRIGHTNESS METHOD 



A number of the Period-Radius relations cited in 
the previous section use the Infrared Surface Bright- 
ness method to determine the angular diameter of 
the Cepheids. However, if mass loss is important 
and causes infrared excess, the the observed infrared 
fluxes will be brighter than the star alone. In this 
case, mass loss will distort the angular diameters pre- 



dicted by near- infrared surface brightness relations from 
iFouque fc GierenI (|1997D , who found 



log 61 = 0.5474 - 


- 0.2^0 + 0.262(y - 


-^)o, 


(4) 


log 61 = 0.5474 - 


-0.2ifo + 0.220(J- 


-K)o. 


(5) 



The IRSB technique is calibrated with individual 
Cepheids, and if these Cepheids have mass loss causing 
an infrared excess, this will distort the predictions of the 
methods. Specifically, four Cepheids used in the calibra- 
tion, V Gen, U Sgr, RZ Vel, and U Car, are predicted to 
have mass-loss rates ra nging from lO"^'^ to 10~^MQ/yr 
(jNeilson fc Lesterll2008f ). These mass-loss rates could in- 
crease the uncertainty of the coefficient multiplying the 
color term, causing a greater scatter in the fitting of the 
calibration sources, and possibly causing the coefficient 
to be overestimated. 

The infrared excess from mass loss will also affect the 
value of the resultant angular diameter in two ways. 
First, it will cause the Cepheid to appear brighter in the 
J-band, and even more so in the X-band. Second, the 



Mass Loss in LMC Cepheids II. 



TABLE 2 
Period-Radius Relations From the Literature 



Reference 



Slope 



Zero Point 



Gieren et aL (1998) 
^ono_^t_^l^ Q998), Canonical 
^ono_^t_^l^ Q998), Non-Canonical 
Gieren et al. (1999) 
Kervella et al. ('2004') 



Ruoppo ct al. (2004) 
Grocncwcgcn (2007) 



0.750 ± 0.024 
0.666 ± 0.007 
0.653 ± 0.006 
0.680 ±0.017 
0.767 ±0.009 

0.69 ±0.09 
0.686 ± 0.036 



1.075 ± 0.007 
1.192 ± 0.009 
1.183 ± 0.009 
1.146 ±0.025 
1.091 ±0.011 

1.22 ±0.08 
1.134 ±0.034 



brighter J- and K- magnitudes will increase the {V — K)q 
and (J — K)o colors, further increasing the angular diam- 
eters. Therefore, mass loss tends to predict larger angu- 
lar diameters of Cepheids and hence larger radii when the 
Baade-Wesselink method is employed for the {V—K)o re- 
lation. It is unclear how the {J — K)o relations is affected 
because the infrared excess makes the Cepheid brighter 
in J-band which acts to counter the increase of the color. 

We begin by assessing how mass loss affects the uncer- 
tainty of the coefficient of the color terms in the relations. 
To fit the stellar radius to the observations, we have as- 
sumed a distance to each Cepheids, which is equivalent to 
fitting each Cepheid's angular diameter. The predicted 
radii and angular diameters derived from the mass-loss 
model differ by less than a few percent from the predicted 
radii and angular diameters from fitting only the radius 
to the observations. This difference between the radii 
has a negligible effect on the comparison of the angular 
diameters from the IRSB technique. 

The mass loss also affects the predicted angular di- 
ameter via the J- and K-hand observations. This in- 
fluence is included when the angular diameters of the 
Cepheids are computed from the observed fluxes from 
the Ngcow et al. (2009) data set, and from the predicted 
stellar fluxes found by fltting the mass-loss model, using 
Equations H] and El Figure [9] compares the angular di- 
ameters using the IRSB technique using observed fluxes 
and the angular diameters derived using the predicted 
fluxes with the best-fit angular diameter from the mass- 
loss model, which is the best-fit radius divided by the 
assumed distance to the LMC. The angular diameters 
from the observed fiuxes exhibit significantly more scat- 
ter than the angular diameters from the predicted fluxes. 
There is reasonable agreement between the angular di- 
ameters from the predicted fluxes from Equation 2] and 
[5l although the comparison from Equation [5] has a larger 
range of differences. 

In the right panel of Figure [9l we show the relative 
difference between the angular diameters predicted from 
the observed and predicted V-, J-, and iiT-band fluxes 
and the best-fit angular diameters as a function of pul- 
sation period. The solid lines represent the mean rela- 
tive difference and the dotted lines are the best-fit lin- 
ear relation. The average difference between the angular 
diameter found using the observed {V — K)o and the 
best-fit angular diameter is 3.35%, while the predicted 
stellar {V — K)q gives a difference of 2.26%. Doing the 
comparison of the best-fit angular diameter found with 
the observed (J — K)q, the relative difference is —0.42%, 
while using the predicted {J — K)q gives a relative differ- 
ence of —0.99%. The lower average uncertainty for the 



(J ~K)o comparison is likely due to the mass loss having 
a similar effect on both the J- and K-ha,nd magnitudes, 
partially canceling out the effect of the color term. The 
X-band magnitudes will also be smaller when infrared 
excess is not removed, thereby predicting a smaller angu- 
lar diameter, which explains why the average difference 
is negative. The relative difference between the angular 
diameters is a measure of the uncertainty of the IRSB cal- 
ibration technique, which suggests that the IRSB tech- 
nique overestimates the radius of a Cepheid by a few 
percent because of the correlation present in the coeffi- 
cient multiplying the color. The differences between the 
angular diameters derived from the observed fluxes and 
the best-fit angular diameters measure the uncertainty of 
both the calibration and the uncertainty due to infrared 
excess. The average uncertainty due to infrared excess is 
about 3 to 4.5% on average, although in rare cases the 
uncertainty may be up to 10 to 20%. 

Figure |9] displays a slight period dependence. One pos- 
sible cause of this dependence is if the stellar color is 
not solely related to the effective temperature. There is 
also a deviation from linearity for small periods. This 
deviation could be due to the interaction of the photo- 
spher e with the hydrogen ionization front ([Simon et al.l 
119931) , which would make the surface brightness relation 
(Equations |4] and [5]) nonlinear functions of color. Be- 
cause the ionization front is closer to the photosphere for 
hotter Cepheids, this interaction has a greater effect for 
shorter-period Cepheids, which alters the effective tem- 
perature from that of a non-pulsation star of the same 
mass, luminosity and radius. For longer-period Cepheids 
the interaction is diminished because the ionization front 
is deeper; the effective temperature and color are more 
similar to a non-pulsating star. Because the calibration 
Cepheids used for the IRSB method may be preferen- 
tially shorter period stars, this could lead to systematic 
differences when applied to longer-period Cepheids. This 
period dependence is similar to the d ependence of the dis - 
tance modulus of the LMC found by lGieren et al.l ()2005f ) 
using a constant value of the projection f actor. In or- 
der to obtain a constant distance modulus, iGieren et al.l 
([2005") derived a period-projectio n factor relation that is 
different from that determined bv lNardetto et ahl ([2007f ) 
from spectroscopy. It is possible that mass loss and IR 
excess contribute to the difference. 

8. PREDICTIONS OF THE LINEAR PERIOD-LUMINOSITY 
RELATIONS 

We have shown that mass loss leads to infrared ex- 
cess that changes the J- and iiT-band fluxes of Cepheids, 
which, in turn, affect the angular diameters determined 
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Fig. 9. — (Left) Comparison of the angular diameters of the Cepheids using the best-fit radius and a distance modulus of 18.5 and the 
angular diameters found by using the IR surface brightness technique of llbugue fc GierenI II1997I ), Equations [4] and [S] using the predicted 
stellar fluxes and the observed fluxes from the average of the two SAGE epochs with a fit x^ < 1-25. The dashed lines represent where the 
predicted angular diameter equals the angular diameter from the IR surface brightness method. (Right) The fractional difference between 
the IR surface brightness technique and the best-fit radius as a function of pulsation period. The solid horizontal lines represent the average 
fractional difference of the sample and the dotted lines (in the plots using the predicted stellar fiuxes) are best-fit linear relations. 



from the IRSB method. The infrared excess may also 
affect the infrared P-L relations. 

In investigating the infrared P-L relations, we continue 
to work with the VIJHK data from two SAGE epochs 
as well as the IRAC data. Because the infrared excess 
due to mass loss will not contribute flux to the V- and 
/-bands, the slopes of the predicted V- and /-band P-L 
relations should be the same as the observed relations. 
However, the structure of the P-L relations will be influ- 
enced by the value of x^ chosen for the threshold; here 
we take x^ < 1.25 as our limit. The zero-points may 
vary because the V- and /-band fluxes depart from per- 
fect blackbodies, but the uncertainty is about 0.1 mag, 
shown in our previous article. The predicted linear re- 
lations are computed using a 3cr-clipping algorithm for 
the SAGE data from each epoch separately as well as 
the average of the two epochs; these are shown in Ta- 
ble 3. There is not a large difference between the pre- 
dicted slopes and zero-points found by fitting the radius 
only and those determined from the mass-loss model, but 
these differences may be attributed to the differences be- 
tween the best-fit radius from each model as seen in the 
P-R relations shown in Table [TJ 

To test if the signature of the infrared excess in the 
P-L relation is impacted by the SAGE detection limit, 
we have removed the 8.0 /im data and recalculated the 
relations. Consistent with what we have found in previ- 
ous sections, the predicted radii and mass- loss rates are 



not significantly changed, and the slopes of the IRAC 
P-L relations agree with the results in Table |3] to within 
the uncertainty. This suggests that the observed brighter 
fluxes for Cepheids with periods logP < 0.7 are real and 
not solely caused by being near the detection limit of the 
survey. 

About one-half of the Cepheids in the sample are pre- 
dicted to have mass-loss rates (> 10"-^** MQ/yr). To 
test the impact of these Cepheids on the IRAC P-L re- 
lation, they have been removed from the sample and 
the P-L relations have been recomputed for the aver- 
age of the two SAGE epochs; the results are Hsted in 
Table [H wh ere no period cut h as been applied to the 
data, as in iNgeow et aP ()2009n . T he IRAC P-L re- 
lations are si milar to the resu lts of INgeow fc Kanbuii 
(2008) and N geow et all (|2009f l. with the slopes having 
the same wavelength dependence, becoming shallower 
with increasing wavelength. The similarity is due to the 
influence of the uncertainty of the IR fluxes caused by 
the unknown pulsation phase in the fitting process in 
the mass-loss model, for which a small change of IR flux 
might produce mass-loss rates orders of magnitude dif- 
ferent. The data are likely still being influenced by IR 
excess. 



iFreedman eFall (|2008[ ) and iMadore et cdl (|2009[ ) pre- 
sented IR P-L relations for LMC Cepheids using 
a selected data set that is very different than th e 
IR P-L relations found by INgeow fc Kanbud (|2008f ). 
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TABLE 3 

Best Fit Parameters for Predicted Linear Period-Luminosity Relations 





Radius Model 








Mass- Loss Model 






Band 


Slope 


Zero Point 


Dispersion 


N 


Band 


Slope 


Zero Point 


Dispersion 


N 


Epoch 1 


V 


-2.745 ± 0.039 


17.235 ±0.028 


0.204 


697 


V 


-2.765 ± 0.039 


17.252 ± 0.028 


0.204 


707 


I 


-2.932 ±0.028 


16.420 ± 0.020 


0.148 


698 


I 


-2.947 ±0.028 


16.433 ± 0.020 


0.148 


707 


J 


-3.081 ±0.021 


16.298 ±0.015 


0.112 


701 


J 


-3.099 ± 0.021 


16.313 ±0.015 


0.109 


706 


H 


-3.135 ±0.019 


16.072 ±0.014 


0.100 


702 


H 


-3.147 ±0.019 


16.084 ± 0.013 


0.099 


709 


K 


-3.181 ±0.017 


15.993 ±0.013 


0.091 


703 


K 


-3.193 ±0.017 


16.005 ± 0.012 


0.089 


710 


3.6 


-3.240 ±0.015 


15.967 ±0.011 


0.080 


698 


3.6 


-3.250 ±0.015 


15.978 ±0.011 


0.079 


705 


4.5 


-3.253 ±0.015 


15.894 ±0.011 


0.079 


698 


4.5 


-3.266 ±0.015 


15.906 ±0.011 


0.078 


706 


5.8 


-3.265 ±0.015 


15.896 ±0.011 


0.077 


698 


5.8 


-3.279 ± 0.016 


15.910 ±0.011 


0.077 


707 


8.0 


-3.275 ±0.015 


15.899 ±0.011 


0.076 


697 


8.0 


-3.288 ± 0.014 


15.913 ±0.010 


0.076 


706 


Epoch 2 


V 


-2.778 ±0.037 


17.265 ±0.027 


0.200 


704 


V 


-2.776 ± 0.036 


17.267 ±0.026 


0.197 


713 


I 


-2.958 ±0.027 


16.440 ±0.019 


0.144 


704 


I 


-2.964 ±0.026 


16.449 ± 0.019 


0.144 


716 


J 


-3.102 ±0.020 


16.314 ±0.014 


0.109 


708 


J 


-3.105 ± 0.020 


16.320 ±0.014 


0.108 


718 


H 


-3.148 ±0.018 


16.083 ±0.013 


0.099 


712 


H 


-3.151 ±0.018 


16.088 ± 0.013 


0.985 


722 


K 


-3.198 ±0.016 


16.005 ±0.012 


0.091 


714 


K 


-3.200 ±0.016 


16.009 ± 0.012 


0.090 


724 


3.6 


-3.251 ±0.015 


15.975 ±0.011 


0.080 


708 


3.6 


-3.251 ±0.015 


15.977 ±0.010 


0.079 


715 


4.5 


-3.260 ±0.014 


15.899 ±0.011 


0.080 


706 


4.5 


-3.261 ±0.014 


15.902 ± 0.010 


0.078 


715 


5.8 


-3.271 ±0.014 


15.900 ±0.010 


0.077 


706 


5.8 


-3.272 ± 0.014 


15.904 ±0.010 


0.077 


715 


8.0 


-3.282 ±0.014 


15.904 ±0.010 


0.076 


706 


8.0 


-3.283 ± 0.014 


15.907 ±0.010 


0.076 


715 


Average of Two Epochs 


V 


-2.720 ±0.034 


17.218 ±0.024 


0.198 


821 


V 


-2.755 ± 0.034 


17.225 ± 0.024 


0.199 


826 


I 


-2.935 ±0.025 


16.422 ±0.018 


0.143 


821 


I 


-2.952 ± 0.025 


16.441 ± 0.018 


0.145 


828 


J 


-3.075 ±0.018 


16.294 ±0.013 


0.108 


825 


J 


-3.099 ± 0.019 


16.318 ±0.013 


0.111 


835 


H 


-3.129 ±0.017 


16.069 ±0.012 


0.099 


831 


H 


-3.144 ±0.017 


16.084 ±0.012 


0.098 


835 


K 


-3.175 ±0.015 


15.990 ±0.011 


0.089 


832 


K 


-3.189 ±0.015 


16.003 ±0.011 


0.089 


836 


3.6 


-3.233 ±0.014 


15.963 ±0.010 


0.079 


826 


3.6 


-3.245 ±0.013 


15.976 ±0.010 


0.078 


828 


4.5 


-3.242 ±0.013 


15.887 ±0.010 


0.077 


825 


4.5 


-3.253 ±0.013 


15.898 ± 0.009 


0.076 


826 


5.8 


-3.256 ±0.013 


15.891 ±0.009 


0.076 


824 


5.8 


-3.269 ± 0.013 


15.904 ± 0.009 


0.076 


829 


8.0 


-3.265 ±0.013 


15.894 ±0.009 


0.076 


825 


8.0 


-3.280 ±0.013 


15.908 ± 0.009 


0.075 


829 



iNgeow et all (I2009D and pre dicted by iNeilson etHI 



()2C)09( ). iFreedman et all (pOOJ ) used SAGE epoch 1 ob 



servations and IRAC data for a sample of Cepheids from 
IPersson et al. (2004) to derive P-L re lations with slopes 
ranging from about -3.30 to -3.44. iMadore et al.l (|2OO90 
used the same Cepheid sample and the average of SAGE 
epoch 1 and 2 observations to derive very similar slopes 
ranging fro m —3.35 to —3 . 49. A s shown in Figure [TOl the 
slopes from lMadore et all ([2009') are significantly steeper 
th an those found her e or by Ngeow & Kanbur (200j) 
or INgeow et al.l (|2009[ ). There are tw o reasons why the 
observed P -L relations differ. First, INgeow fc Kanburl 
((2008) and INgeow et all ()2009D used an iterative fitting 
routine designed to minimize the effect of Gepheids with 
fluxes that are very different from the majority of the 
sample, possibly due to effects such as blending or false 
identification. The method retains all the Gepheids but 
changes their weighting in fitting the P-L relation, thus 
affecting the structure of the IR P-L relations; the slopes 
are shallower. The second reason for the difference is 
th at the period range of t he samples. The Gepheids used 
bv lFreedman et al.l ()2008D and Mado re et all (|2009!) have 
periods ranging from about 6 days to 50 days, with 83% 
of Gepheids having periods greater than 10 day s. The 
sample of Gepheids used by iNgeow et al.l (|2009l ) and in 
this paper are mostly shorter period. To test the sen- 
sitivity of the slope to the periods, we computed the 



IRAG P-L relations for our dataset using different min- 
imum periods, after removing those Gepheids with an- 
gular separation > 0.5" between the surveys and also 
removing those Gepheids that are consistent with mass 
loss. We find that when the shorter period Gepheids are 
removed from the sample, the slopes of the IRAG P-L re- 
lations steepen. For a minimum period of logP = 1.05, 
the slopes of the IRAG P-L relations, computed using a 
3(7 iterative method, are c onsistent with the results of 
[fteedman et al. (200l) and IMadore et al.l ([2009); those 
P-L relations are also listed in Table HI 

9. CONSISTENCY OF THE DIFFERENT DATA SETS 

It is important to note that the sample of Gepheids 
used by Freed man et al.l (|2008| ) and IMadore et al.l (|2009f ) 
are consistent with the sample used her e as concerns the 
hypot hesis of Gepheid ma s s loss. In the IFreedman et al.l 
(|2008D and IMadore et all (|2009[ ) samples there are no 
Cepheids with significant IR excesses. This is also true of 
the sample used here. Of the 200 Gepheids with P > 10 
days having x^ < 1-25 in the average SAGE epoch that 
we analyzed for the statistical significance of mass loss, 
we find only two Gepheids that are consistent with mass 
loss. From t he consistency o f our s ample and the sam- 
le us ed by IFreedman et all (|2008D and IMadore et al.l 
2009f) . and from the inherent difficulty of determining 
mass loss from the data of long period Cepheids, we sug- 
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The difference between the IR P-L relations of IJNgeow et al.l 120091) and IMadore et al.l 120091 ) and the IR P-L relations in this work. The 
horizontal line represents zero difference. 



gest that these samples are aU affected by mass loss to 
the same degree. Taking period cuts of the predicted IR 
fluxes and computing the best-fit IR P-L relations pre- 
dicts slopes that become shallower with increasing pe- 
riod. The observed steeper P-L relations defined by the 
sample of lon ger period Cepheids is c onsistent with the 
predictions of iNeilson fc Lesten ()2009t ) that mass loss for 
these Cepheids is driven primarily by radiation, with- 
out major contributions from shocks and pulsation, with 
the mass-loss rates ranging from 10^^ — lO^^MQ/yr. 
Radiative-driven mass loss is period dependent, which 
leads to an infrared excess compared to the stellar lumi- 
nosity for Cepheids with P > 10 days. Therefore mass 
loss might cause the steeper IR P-L relations. 

Another consistency check is the corresponde nce be- 
tween the P-R and the P-L relations. Freedm an et al.l 
(|2008h were able to derive the slope of the P-L relation 
in the limit of long wavelengths by using the result from 
iGieren et al.l ()1999f l that the slope of the P-R relation is 
about 0.68. Because the long wavelength flux is primar- 
ily dependent on the Cepheid's radius, they ware able to 
flnd that the slope of the P-L relations is about —3.4. 
In this work, we find that the slop e of th e P-R relation 
is the same as that of IGieren et al.l (|1999f ) but the slope 
of the predicted IR P-L relations range from -3.3 to -3.2 
significantly different from that asymptotic limit. We in- 
vestigate the apparent contradiction by re-deriving the 
asymptotic limit to the slope of the P-L relation and 
plotting the predicted and observed 8.0 ^m fluxes as a 
function of the predicted radius in Figure 1111 In Equa- 
tion [1] we had that L^ ex E?Bv{Tes): which means the 



flux is fi, oc R^B^{Tcs)/cP. The distance d is that of the 
LMC and is assumed to be the same for all Cepheids in 
our sample. Converting the flux to apparent magnitude 
yields 



mA = -51ogi?- 2.51ogP,(reff) -f Ci, 



(6) 



where C'l is a constant and contains the distance. If the 
slope of the linear best fit to the magnitude versus the 
log i? is — 5 then there is a serious problem with our anal- 
ysis as it suggests that the Blackbody function does not 
contribute to the brightness in the infrared. However, 
the best-fit relations have slopes of —4.70 and —4.12 for 
the predicted and observed fluxes respectively. This sug- 
gests that the radius of the Cepheid is related to the 
Blackbody radiation via the effective temperature and 
pulsation period. The difference between the slopes cor- 
responds to th e slopes of th e resp ective P-L relations: 
-3.49 from the IMadore et al.l ()2009f ) sample, -3.28 from 
our predicted fluxes, and —2.85 from the observed fluxes. 
As expected the flux-radius relations do not have a slope 
of —5 because the effective temperature affects the fit 
of both the predicted and observed fluxes, even at long 
wavelengths via its relation to the radius. Using the ob- 
served fluxes, the apparent IR excesses cause the slope 
to become even shallower. We show this by replacing 
the blackbody function at infrared wavelengths by the 
Rayleigh- Jeans approximation, B^(Tcg) ex Toff, giving 

mA = -51ogi?-2.51ogTeff + C2. (7) 

iBeaulieu et al.l (|2001f ) have shown that the effective tem- 
perature is a function of the gravity and {V — /)o color. 
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TABLE 4 
Best Fit Parameters for Observed IRAC Linear Period-Luminosity Relations 



Band 


Slope 


Zero Point 


Dispersion 


N 




Observed IR Fluxes With Mass-Losing Cepheids 


Removed 


3.6 


-3.224 ±0.015 


15.984 ±0.011 


0.089 


761 


4.5 


-3.190 ±0.015 


15.945 ±0.011 


0.089 


763 


5.8 


-3.203 ±0.020 


15.938 ±0.015 


0.120 


763 


8.0 


-3.080 ±0.033 


15.818 ±0.028 


0.134 


429 


Observed IR Fluxes With logP > 1.05 Mass-Losing Cepheids Removed 


3.6 


-3.239 ±0.160 


16033 ±0.191 


0.135 


62 


4.5 


-3.068 ±0.172 


15.849 ± 0.206 


0.146 


63 


5.8 


-3.168 ±0.166 


15.931 ±0.198 


0.140 


63 


8.0 


-3.256 ±0.175 


16.016 ±0.175 


0.148 


63 


Observed IR Fluxes 


3.6 


-3.266 ±0.010 


15.991 ±0.007 


0.100 


1700 


4.5 


-3.210 ±0.010 


15.938 ±0.007 


0.102 


1715 


5.8 


-3.072 ±0.017 


15.775 ±0.012 


0.158 


1521 


8.0 


-2.854 ±0.031 


15.500 ±0.026 


0.221 


815 


Observed IR Fluxes for logP > 1.05 


3.6 


-3.347 ± 0.068 


16.112 ±0.087 


0.150 


118 


4.5 


-3.320 ±0.067 


16.093 ±0.087 


0.150 


118 


5.8 


-3.395 ±0.069 


16.42 ± 0.089 


0.153 


119 


8.0 


-3.517 ±0.074 


16.264 ±0.096 


0.165 


119 
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Fig. 11. — The (Left) predicted and (Right) observed 8.0 /im flux of the LMC Cepheids as a function of the predicted radii of the 
Cepheids. The solid line is the linear least squares best-flt to the data, while the dotted line is the linear least squares best-fit writh the 
slope forced to be —5, which is the expected slope if the temperature of the star does not effect the flux. 



logTeff = 3.92 -H 0.0055 logg - 0.2487(V^ - /)o. The 
color term can be replaced by the period via a Period- 
Color (P-C) relation and the gravity term is \ogg = 
2.62 — 1. 21 log P. For i nstan ce, using the P-C rela- 
tion from iSandage et"aLl ()2004[ ) for the LMC Cepheids, 
{V -I)q = 0.256 log P 4- 0.444, gives 



mx 



-51ogi? + 0.1751ogP-|-C3 



(8) 



If we use the P-R relation to replace the period term, we 
can derive the limiting slope to the flux-radius relation. 



mx 



-4.74 log i? + Constant. 



(9) 



Or we can use the P-R relation to replace the radius term 
to find the limiting slope of the P-L relation. 



mx 



-3.231ogP + Constant. 



(10) 



The slope of Equation |9] agrees with the best-fit slope 
of the flux-radius relation from Figure llll suggesting 
that indeed the effective temperature affects the slope 
of the P-L relation at long wavelengths. Therefore 
the asymptotic limit of the slope of the P-L relation 
at long wavelengths is approximately —3.23. However, 
this does not neces sarily mean that the P-L relations 
of lFreedmanet all (|2008| ) and iMadore et all (|2009| ) are 
incorrect. If the slope of the P-R relation is 0.71, the 
asymptotic limit of the slope of the P-L relation is —3.4, 
still consistent wit h ob served P-R rela t ions found by 
iGieren et all (|1998D and lLanev fc Stobid (|1995f ). 

10. PREDICTIONS OF THE NON-LINEAR 
PERIOD-LUMINOSITY RELATIONS 



iNgeow fc Ka^Aml poM) and iNgeow et al] (|2009l ) 
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found the surprising result that at 5.8 and at 8.0 fim the 
sfopes of the P-L relations become steeper for Cepheids 
with P > 10 days. We want to test if this could be 
caused by mass loss that makes the predicted P-L rela- 
tions non-linear. To do this, we have used two linear P-L 
relations, one for periods less than 10 days and the other 
for longer periods. These are of the form 



m\ 



alogP + b 
clogP + d 



log P < 1 
logP> 1 



(11) 



We compute the best-fit non-linear relations employing 
the same Scr-iterative routine used before. The F-test is 
used to check if the non-linear relations are statistically a 
better fit to the data and the relations are listed in Table 
[H We find that all P-L relations at all wavelengths for 
each epoch have P > 3, meaning they are consistent with 
being non-linear. Mass loss affects the structure of the 
P-L relations and the potential value of F computed for 
the IR P-L relations. This indicates that the presence 
of mass loss acts to hide the non-linear nature of the IR 
P-L relations, steepening the slopes for logP > 1 of the 
observed 5.8 and 8.0 /im P-L relations. All the predicted 
P-L relations have slopes for log P > 1 that are shallower 
than the slopes for logP < 1 for the same wavelength 
for second and average epoch fluxes. The first epoch 
fiuxes have non-linear P-L relations where the slopes for 
the long-period Cepheids is approximately the slope of 
the P-L relation for the short-period Cepheids within the 
uncertainty. Because the errors of the slope of the long- 
period relations is large then the structure of the first 
epoch non-linear P-L relations are still consistent with 
the structure of the non- linear P-L relations for the other 
epochs. This also shows why the slopes of the predicted 
IR linear P-L relations increases from —3.3 to —3.2 when 
the short-period Cepheids are removed from the sample, 
while the opposite effect seen for the observed Cepheid 
fiuxes. 

In section 9 we based our analysis on the asymptotic 
limit of the slope of the linear IR P-L relations. Here we 
want to use this limiting behavior to test if it is plausible 
that the IR P-L relations are actually non-linear. In 
section 9 we used Equation [T] and the Rayleigh- Jeans 
approximation for the blackbody function, which enabled 
use to write the effective temperature as a function of 
the {V — I)o color and then as a function of period via 
the Period-Color relation. Here we replace the linear 
Period-Color relation with the non -linear Period-Color 
relation from iSandage et al.l (|2004[ ). which has a slope 
of 0.160 for logP < 1 and a slope of 0.315 for logP > 
1. Doing this readily shows that the IR P-L relation is 
non-linear at long wavelengths. Furthermore, the P-L 
relation is the sum of a linear P-R relation and a non- 
linear Period-Color relation. In that case the limiting 
slope of the non-linear P-L relation is -3.30 for logP < 
1 and the limiting slope is -3.21 for logP > 1. These 
values are consistent with the non-linear IR P-L relations 
we have derived here, but we note that the analysis is 
very sensitive to the slope of the P-R and Period-Color 
relations. An increase of only 0.01 to the slope of the 
P-R relati on causes the P-L relation s to be steeper by 
0.05. Also. lKanbur fc Ngeowl (|2004bf ) found a non-linear 
Period-Color relation for LMC Cepheids with a slope of 
0.152 for logP < 1 and a slope of 0.590 for logP > 1. 



The limiting slopes of the non-linear IR P-L relations is 
then 3.31 and 3.04. 

Our derivation also provides insight into why 
iNgeow et al.l (l2009[ ) found the IR P-L relations that are 
consistent with linearity. At any wavelength the P- 
L relation is the combination of the P-R relation and 
the Period-Color relation, connected via the effective 
temperature. The P-R relation has a slope and zero- 
point that are independent of wavelength, but the effec- 
tive temperature contributes to the luminosity through 
the wavelength-dependent blackbody function. At wave- 
lengths in the Rayleigh- Jean regime, the blackbody func- 
tion is a linear function of the effective temperature, 
while at optical wavelengths the greater dependence of 
the full Planck function on effective temperature leads 
to a more dramatic change in luminosity, i.e. the slope 
of non- linear Period-Color relation has a much more im- 
portant contribution to the structure P-L relation. At 
longer wavelengths, the non-linear Period-Color relation 
plays a much reduced role, and the change of the slope 
with period is minimum. Thus, to be able to determine 
if the IR P-L relations are non-linear the scatter of the 
fit of observed fiuxes of Cepheids must be significantly 
smaller and have any contamination of IR excess from 
mass loss removed. 

Our predicted linear P -L relations are similar t o the 
observed relations from INgeow fc Kanbuil (|2008D and 
Ngeow et al. (2009), suggesting that mass loss does not 
contribute much to the IR observations. However, the 
non-linear IR relations are significantly different. For 
logP > 1 the relations have predicted slopes that are 
shallower than the observed relations, while for log P < 1 
the predicted slopes are steeper than the slopes of the 
observed relations. These results suggest that mass loss 
may play two roles, as can be seen in Figure 1121 which 
plots the ratio of the fiux from the dust shell to the total 
fiux of the Cepheids as a function of the pulsation period 
for each IRAC band. The first role of mass loss is to 
increase the dispersion of the P-L relations at all period 
cuts, and the second role is to alter the slope of the linear 
and non-linear P-L relations. 

The dispersion of the IR P-L relations for LMC 
Cepheids is caused by several factors: the width of the 
Cepheid instability strip, the tilt and depth of the LMC 
and possibly mass loss. Considering the geometry and 
the metallicity distribution of the LMC, iMadore et al.l 
(2009) derived the minimum observable dispersion of the 
linear IR P-L relations to be about 0.08 mag. Our anal- 
ysis here finds the dispersions of the linear IR relations 
to range from 0.079 mag at 3.6 /xm to 0.074 mag at 8.0 
/xm. These dispersions are consistent with the minimum 
dispersion found by Madore et al. (2009) , which suggests 
that mass loss contributes only 0.001 to 0.006 mag to the 
dispersion of the IR P-L relations. 

The behavior of mass loss in LMC Cepheids affects 
the structure of the IR P-L relations. As shown in 
Figure fT2l for shorter-period Cepheids (logP < 1) the 
shell flux is a larger fraction of the total flux, while for 
longer-periods Cepheids the ratio is approximately con- 
stant within some dispersion. The large number of short- 
period Cepheids in the sample causes the slopes of the 
linear IR P-L relations to be shallower, and as wavelength 
increases, these slopes become even more shallow. This 
behavior explains why the slopes of the IR P-L relations 
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TABLE 5 

Best Fit Parameters for Predicted Non-Linear Period-Luminosity Relations 



Band 



Slope5 



Zero Point c 



Slopes 



Zero Pointy Dispersion N 



Epoch 1 



V 


-2.917 ±0.056 


17.344 ± 0.037 


-2.457 ±0.273 


16.958 ±0.313 


0.202 


707 


7.09 


I 


-3.068 ±0.041 


16.507 ± 0.027 


-2.743 ±0.196 


16.246 ±0.225 


0.146 


707 


8.59 


J 


-3.185 ±0.030 


16.366 ± 0.020 


-2.932 ±0.139 


16.155 ±0.160 


0.108 


708 


7.87 


H 


-3.218 ±0.027 


16.127 ±0.018 


-3.113 ±0.125 


16.073 ±0.143 


0.098 


709 


6.33 


K 


-3.255 ±0.025 


16.042 ±0.016 


-3.162 ±0.113 


15.995 ±0.130 


0.089 


710 


5.99 


3.6 


-3.296 ±0.022 


16.006 ±0.015 


-3.281 ±0.101 


16.033 ±0.116 


0.078 


705 


4.59 


4.5 


-3.310 ±0.021 


15.933 ±0.014 


-3.304 ±0.100 


15.970 ±0.114 


0.078 


706 


4.39 


5.8 


-3.321 ±0.021 


15.936 ±0.014 


-3.325 ±0.099 


15.981 ±0.113 


0.077 


707 


4.13 


8.0 


-3.327 ±0.021 


15.936 ±0.014 


-3.343 ± 0.097 


15.993 ±0.111 


0.076 


706 


3.68 


Epoch 2 


V 


-2.928 ±0.053 


17.360 ± 0.035 


-2.501 ±0.208 


17.005 ±0.241 


0.195 


713 


7.73 


I 


-3.077 ±0.039 


16.519 ±0.26 


-2.767 ±0.152 


16.263 ±0.176 


0.143 


716 


8.09 


J 


-3.190 ±0.029 


16.372 ±0.019 


-2.945 ±0.111 


16.165 ±0.129 


0.107 


718 


8.08 


H 


-3.219 ±0.027 


16.129 ±0.18 


-3.081 ±0.101 


16.032 ±0.116 


0.098 


722 


5.94 


K 


-3.255 ± 0.024 


16.043 ± 0.016 


-3.147 ±0.092 


15.968 ±0.107 


0.090 


724 


4.66 


3.6 


-3.299 ±0.021 


16.006 ±0.014 


-3.214 ±0.081 


15.954 ±0.094 


0.079 


715 


4.60 


4.5 


-3.311 ±0.021 


15.933 ±0.014 


-3.197 ±0.081 


15.848 ± 0.093 


0.077 


715 


5.37 


5.8 


-3.320 ±0.021 


15.933 ±0.014 


-3.213 ±0.080 


15.853 ±0.092 


0.076 


715 


4.85 


8.0 


-3.327 ±0.021 


15.934 ±0.014 


-3.227 ±0.079 


15.861 ±0.091 


0.076 


715 


4.36 


Average of Two Epochs 


V 


-2.882 ±0.050 


17.329 ± 0.033 


-2.429 ±0.210 


16.921 ±0.245 


0.198 


826 


6.56 


I 


-3.045 ± 0.036 


16.498 ± 0.024 


-2.775 ±0.153 


16.273 ±0.177 


0.144 


828 


6.45 


J 


-3.170 ±0.028 


16.361 ±0.018 


-2.921 ±0.114 


16.137 ±0.133 


0.110 


835 


6.72 


H 


-3.200 ±0.025 


16.118 ±0.016 


-3.066 ±0.100 


16.016 ±0.116 


0.098 


835 


4.98 


K 


-3.238 ±0.022 


16.033 ±0.015 


-3.121 ±0.091 


15.944 ±0.106 


0.089 


836 


4.66 


3.6 


-3.284 ±0.020 


15.999 ±0.013 


-3.211 ±0.080 


15.952 ± 0.093 


0.078 


828 


3.61 


4.5 


-3.290 ±0.019 


15.921 ±0.013 


-3.200 ±0.078 


15.852 ±0.091 


0.076 


826 


3.53 


5.8 


-3.309 ±0.019 


15.928 ±0.012 


-3.217 ±0.078 


15.860 ±0.091 


0.076 


829 


4.03 


8.0 


-3.316 ±0.019 


15.930 ±0.012 


-3.233 ± 0.077 


15.869 ±0.090 


0.075 


829 


3.60 



determined by iNgeow et al.l (|2009[ ) become shallower as 
a function of increasing wavelength. Mass loss likewise 
affects the structure of the non-linear IR P-L relations. 
For the short period part of the non-linear IR P-L rela- 
tions, mass loss causes the slope to be shallower, while for 
the long period Cepheids, the slope is steepened. At 3.6 
^m and 4.5 fim, this, along with an increased dispersion, 
makes the P-L relations appear to be consistent with lin- 
earity by causing a smaller computed value of F for the 
F-test. At longer wavelengths this effect is more pro- 
nounced; mass loss causes the slope of the P-L relations 
for logP < 1 to be more shallow than the long-period 
relations, opposite of what has been determined for the 
optical and near-IR relations. 

The predicted linear and non-linear P-L relations, 
as well as the effects of mass loss on their structure, 
agree with the results of our previous work using SAGE 
epochl-data correlated with the OGLE-11 Cepheids. In 
that work the IR P-L relations were found to be shal- 
lower than those found here, the difference being due to 
the use of a stricter x^ cutoff, which we noted has an 
effect on the predicted slope of the P-L relations. Re- 
gardless of this difference, we have shown in both works 
that the IR P-L relations are affected by infrared excess 
caused by the presen ce of d ust shells similar to those that 
have been observed (iKe rvclla et al. 2006; iMerand et al I 
l2006ll2007tlMarengoetral...2009.) . Because the amount of 
circumstellar dust is dependent on the dust-to-gas ratio, 
the infrared excess is metallicity dependent. The IR P-L 
relations, therefore, are metallicity dependent as shown 



in the previous work. 

11. CONCLUSIONS 

The purpose of this work is to use the OGLE-III sample 
of fundamental-mode Cepheids, complemented by data 
from 2MASS and two epochs of SAGE observations from 
INgeow et al.l ([2009) to test for the existence of infrared 
excess due to the formation of dust shells at a large dis- 
tance from the surface of a Cepheid. The dust shells 
are hypothesized to be caused by mass loss from the 
Cepheids. In addition to testing for the existence of IR 
excess, we have investigated how the IR excess affects the 
structure of the P-L relations, whether the P-L relations 
are consistent with being non-linear and the uncertainty 
of the infrared surface brightness technique due to mass 
loss. 

We fit two models to the VIJHK and IRAC observa- 
tions; the first model fits the observed fiuxes of Cepheids 
using just one free parameter, the radius of the Cepheids. 
The second model uses two free parameters, the radius 
and the dust mass-loss rate which is converted to a gas 
mass-loss rate. In both cases, the effective temper ature is 
calculated using the relation from lBeaulieu et al.l ([2001). 
The fits of the two models are compared using the F- 
test, and we found that the sample is consistent with 
mass loss. In addition to having large mass-loss rates, 
these Cepheids tend to have shorter periods, logP < 1. 
This is because long-period Cepheids are cooler and have 
larger IR stellar fluxes, which makes the contribution of 
a circumstellar dust shell a smaller relative contribution 
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Fig. 12. — The predicted ratio of the luminosity of the dust shell surrounding each Cepheid generated by mass loss to the total (shell 
stellar) luminosity of each Cepheid in the four IRAC wavelengths as a function of pulsation period. 



or IR excess. In this work, we also use a stricter x^ cutoff 
than in our previous paper, which leads us to predict V- 
and /-band P-L relations that agree with observations. 
The use of the stricter cutoff removes Cepheids that are 
poor fits due to being at a distance very different from 
the assumed LMC distance modulus of 18.5 and having 
IR fluxes that are significantly different from the mean 
brightness. 

Because mass loss causes an infrared excess, it also af- 
fects the IR surface brightness technique for determining 
the angular diameter of Cepheids by increasing the ob- 
served J- and iiT-band fluxes to be brighter than the stars 
brightness. For the V- and iiT-band relation, the value of 
{V — K)o is larger and the computed angular diameter is 
overestimated. The error of the relation is about 3.5%. 
Using the J- and iiT-band relation, there are offsetting 
errors because it uses the color (J — K)o and the K- 
band brightness, and has an uncertainty of about 0.4%. 
Also, it is found that the calibration of the IRSB method 
has error due to the possibility of mass loss in the cali- 
brating Cepheids. The uncertainty due to calibration is 
about 2.2% and 0.9% for the V, K-hand and J, X-band 
relations, respectively. 

Using our predicted radii, we calculated the Period- 
Radius relations for each epoch of SAGE observations. 
We find that the P-R relations have a slope of 0.690 ± 
0.003 and a zero-point of 1.128 ± 0.002. The slope and 
zero-point agree wi th the slope and zero-point found by 
iGieren et alJ (|1999f ): however the zero-point in our rela- 
tion is notably smaller. 

We predicted the linear P-L relations for the LMC 
Cepheids, finding that the IRAC P-L relations have 
slopes that range from —3.22 to — 3.26. These slopes 
agree with the results of iNgeow et al . (2009) within the 
disper sio n, but not with the results of .Freedman et all 
(|2008f ) or lMadoreet all (|2009D . Because the F-test pre- 
dicts that about half of the Cepheids are consistent with 
no mass loss, we compute the IR P-L relations for the 
sample of Cepheids after removing those that appear to 
have mass loss as indicated by an infrared excess. These 
relations are computed without the period cut-offs sug- 



gested bv INgeow et al.l ()2009[ ). and it is clear the rela- 
tions are not significantly affected by the removal of the 
Cepheids that might have mass loss. This suggests that, 
although the IR excess might still important in the sam- 
ple, it cannot be statistically identified because the un- 
certainty of the observations is so large. The non-linear 
P-L relations are calculated and tested using the F-test, 
and the P-L relations are non-linear at all wavelengths. 
This result differs from the observations where the rela- 
tions are found to be non-linear at wavelengths shorter 
than K-hand. 

We found that the IR excess due to mass loss is strongly 
influenced by the unknown phase of the IRAC fluxes. For 
instance, if the amplitude of light variation for a Cepheid 
in the IR is 0.1 magnitude then the ratio of the maximum 
luminosity to the mean luminosity LmaxI Lmean ~ 1-1 
This would appear to be equivalent to an IR excess from 
a shell Lsheii/Lmean ~ 10%. A comparison to the shell 
fluxes shown in Figure [12] suggests that the unknown 
pulsation phase of the IR fluxes is the most important 
uncertainty in our analysis. Time-series observations 
of Cepheids in the IR with complete phase coverage is 
necessary to better understand how mass loss works in 
Cepheids and how the IR excess affects the IR P-L rela- 
tion and IRSB technique. 

All of these results suggest that mass loss in Cepheids 
has an important effect on the structure of the IR Period- 
Luminosity relations and the Period-Radius relations 
that are predicted from the IR surface brightness tech- 
nique. We have shown that mass loss occurs in half of the 
sample of OGLE-III Cepheids, making it difficult to state 
conclusively that mass loss is an important phenomena 
in all Cepheids. If the uncertainty of the fluxes were less, 
more Cepheids might be consistent with having a stel- 
lar wind. While having a larger sample of Cepheids and 
two epochs of SAGE observations has provided a bet- 
ter testbed for studying mass loss, even more data are 
needed. Fitting mass loss and radii of Cepheids is better 
done with time series observations of the Cepheids to de- 
termine the mean IR flux and decrease the phase uncer- 
tainty. The uncertainty of the IR brightness of Cepheids 
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due to mass loss needs to be characterized il the IR P-L 
relations are to be used as tools for measuring precise dis- 
tances to galaxies and determining the Hubble constant 
to an accuracy of 2-3%. This issue may only be resolved 
by more observations of Cepheids in IR wavelengths. 



HN would like to thank Dr Tom Barnes, and Dr Jes- 
per Storm for enlightening discussions about IRSB tech- 
nique. We would also like to thank the anonymous ref- 
eree for helpful comments. 
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